Spectroscopic signatures of nonequilibrium pairing in atomic Fermi gases 
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We determine the radio-frequency (RF) spectra for non-stationary states of a fermionic condensate 
produced by a rapid switch of the scattering length. The RF spectrum of the nonequilibrium state 
with constant BCS order parameter has two features in contrast to equilibrium where there is a 
single peak. The additional feature reflects the presence of excited pairs in the steady state. In the 
state characterized by periodically oscillating order parameter RF-absorption spectrum contains two 
sequences of peaks spaced by the frequency of oscillations. Satellite peaks appear due to a process 
where an RF photon in addition to breaking a pair emits/absorbs oscillation quanta. 
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Cooper pairing in ultra-cold Fermi gases has been a 
major focus of research in the past few years. Remark- 
able experimental techniques such as sweeps across the 
Feshbach resonance [1 J^J , generation of collective modes 
[H and vortex lattices J and radio-frequency (RF) spec- 
troscopy @, 0] have been developed to probe the paired 
state. While it was crucial to establish for cold gases well- 
known signatures of fermionic pairing, of a key interest 
are regimes not easily accessible in superconductors, e.g. 
strong interactions in the vicinity of the Feshbach reso- 
nance and highly imbalanced mixtures. 

One of the most interesting possibilities is to access 
the non-adiabatic coherent dynamics of fermionic con- 
densates Driven out of equilibrium by a sudden 
change of the pairing strength on the BCS side of the 
Feshbach resonance, these systems acquire steady states 
with properties strikingly different from equilibrium ones. 
Three distinct non-stationary states have been predicted 
- a state where amplitude of the BCS order parameter 
A(t) oscillates periodically 0, [3], a state with a con- 
stant but reduced gap, and a gapless superfluid state 0- 
[l6j . Realization of a particular steady state is determined 
by the magnitude of change of the pairing strength. Most 
previous studies concentrated on the time evolution of 
the order parameter, while direct experimental manifes- 
tations of the non-adiabatic dynamics have not been suf- 
ficiently explored. The purpose of the present paper is 
to address this issue. 

Amongst existing experimental techniques the RF 
spectroscopy appears to have the greatest potential for 
distinguishing different dynamical states from equilib- 
rium phases. This motivates us to study spectroscopic 
signatures of the dynamics of fermionic condensates. Our 
main findings are as follows. In contrast to the BCS 
ground state spectrum which has a single peak at a fre- 
quency determined by the equilibrium gap, the RF spec- 
trum of a nonequilibrium state with constant but finite 
A s displays two distinct peaks. The second peak reflects 
the fact that this nonequilibrium state is a superposition 
of an infinite number of excited stationary states of the 
condensate. Excited states contain a mixture of ground 



state pairs and excited pairs - two-particle excitations of 
the condensate that conserve the total number of par- 
ticles and Cooper pairs (see 17, and the discussion 
below) . 

The "ordinary" peak present already in the ground 
state is due to a process whereby a photon breaks a 
ground state pair, while in the process responsible for 
the second peak it breaks an excited pair, (Fig. [2(a)). It 
is interesting to note that in electronic superconductors 
excited pairs carry no charge or spin and are therefore 
difficult to detect. 
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FIG. 1: (a) Non-stationary state of a fermionic condensate 
with time-independent order parameter A 3 contains a mix- 
ture of ground state and excited Cooper pairs of atoms |1) 
and 1 2} with energies ±\f e 2 + Aj. An RF photon can break 
either type of pair and transfer one of the atoms to state |3). 
(b) In the steady state where the order parameter A(f.) os- 
cillates with frequency Q, the photon can break a pair and 
emit/absorb several oscillation quanta. 



For a steady state with periodically oscillating order 
parameter, we show that each of the peaks described 
above acquires equidistant satellite peaks, i.e. there are 
two series of equally spaced peaks in this state. The spac- 
ing between peaks in each series is equal to the frequency 
of oscillation ft. Satellite peaks appear because a photon 
can gain optimal energy for breaking a ground/excited 



2 



pair by emitting or absorbing several "deltons" - oscilla- 
tion quanta of energy f2, (Fig. HJb)). 

In an atomic Fermi gas the pairing occurs between 
atoms in two hyperfine states |1) and |2). The frequency 
of external RF radiation can be tuned to induce transi- 
tions between one of these states, say |2), to the third 
atomic state |3). The RF spectrum corresponds to the 
rate of loss of atoms from |2), i.e. /(w r f) = — dN^/dt, 
measured as a function of the radiation frequency w r f . In 
the normal state of atoms |1) and |2) the quantity /(w r f) 
has a sharp peak at w r f = u> a , the frequency of atomic 
transition between |2) and |3). In the paired ground state 
the peak shifts to a larger frequency since now an addi- 
tional energy is required to break pairs dEi. 

We start with the Hamiltonian H = H^f s + H 3 + 



H^ M , where 



H 



BCS 
12 



3,a=l,2 



£ j c ja C O a 



Vp 



i,3 



is the BCS Hamiltonian describing pairing between states 
|1) and |2), Cj a (a = 1,2) annihilate atoms in states |1) 
and 1 2), £j are single-particle energy levels relative to the 
Fermi level of atoms |1) and |2), X(t) and vp are the 
dimensionless coupling and the density of states at the 
Fermi level, H% = JT Sjdjdj, where dj annihilate atoms 
in states |3), represents non-interacting atoms in states 
|3), and 
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accounts for the interaction of atoms |2) and |3) 
with the RF radiation field[19| in the rotating wave 
approximation [2^|. Here to = c<J r f — u) a is the detuning 
frequency. Since the size of the trap is much smaller 
than the photon wavelength, one can take the tunnelling 
matrix to be diagonal, Tji — TSji. 

We assume that the pairing strength has been switched 
from Xi to A/ and the RF radiation is turned on after the 
condensate has reached one of the steady states described 
above. The magnitude of the change in pairing strength 
is denoted by the parameter (3: 

f3 = \r l -\f. 

Our task is to evaluate the current (J) = —d(N 2 ) /dt. The 
wave function of the condensate in the steady state with- 
out the RF field is of the BCS form |*(f)) = EL/M*) + 
u j(t)CjiCj2\ |0) [SI- Treating the tunnelling Hamiltonian, 
Ht in eq. ([2]), as a perturbation, we obtain the current 
out of state |2) to the lowest nonvanishing order in Tji 
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FIG. 2: RF spectra © for a non- stationary state with a con- 
stant order parameter A s 7^ produced by an abrupt change 
in the pairing strength, Ai — > A/ for Sfi = —0.75 A/; (a) The 
spectral weight for u < where the peak is due to processes 
where the photon breaks an excited pair. This peak is absent 
in the ground state, = 0; (b) The peak at ui > is due to 
processes where an RF photon breaks a ground state Cooper 
pair. A similar peak is present in the paired ground state. A / 
is the equilibrium gap for the final coupling A/ 



where Uj(u>) are Fourier components of Uj(t) and 5 fx = 
M3 — A*2 is a difference between the corresponding chem- 
ical potentials for atoms in states |3) and |1), |2). In 
recent experiments all states |3) were initially unpop- 
ulated d, which suggests that we set Sfj, ~ — Ep. 
However, our model based on truncated BCS Hamilto- 
nian (fT]) becomes invalid for that case since the so-called 
" off-diagonal" interaction terms between atoms in states 
|1) and 1 2) can not be discarded Q. To circumvent this 
problem in what follows we assume \Sfi\ -C Ep. 

Consider first the steady state with a constant order 
parameter A s that is realized for — 7r/2 < (3 < ir/2. The 
steady state wave function has been determined exactly 
in Ref.0 
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are the Bogoliubov amplitudes in the BCS ground state 
with gap A s . The distribution function cos 2 [0(ej/2)], 
(Fig. [3]) , and A s are known exactly in terms of the ini- 
tial and final equilibrium BCS gaps Aj and A f [l6[ . The 
first term in eq. ^ is the wave function of a ground state 
pair of energy — £j . The second term is the wave function 



of an excited pair with energy £j 17|. Excited pairs are 
excitations of the condensate and should be contrasted 
to the single-particle excitations, which are created out- 
side of the condensate. When the BCS wave function is 



projected onto the subspace of fixed particle number [2l| . 
excited pairs conserve the total number of paired atoms, 
while quasiparticle excitations break Cooper pairs. In 
the Anderson pseudospin representation[2lj, excited and 
ground state pairs correspond to a pseudospin respec- 
tively aligned parallel or antiparallel to its effective mag- 
netic field. In this case 9j is the angle between the pseu- 
dospin and the field. 

Using eqs. (|3l4p . we derive the rate of loss in state |2) 
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FIG. 3: The probability n(e) = cos 2 [#(e)/2] of having an ex- 
cited pair (see the text below eq. (Q) at energy e in all steady 
states produced by a switch of the BCS coupling constant 
\i — » A/ Plots for three values of — 1/Xi — 1/A/ 

are shown. The presence of excited pairs leads to additional 
peaks in RF spectra shown in Figs. 121 II In the ground state 
n(e) = 0. 
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where 



= ^Sfi 2 + A 2 ± 5/j, and To = {uj 2 - A 2 )/2w. 
The first term represents the contribution of ground state 
pairs, (Fig. HJb)), corresponding to a process where a 
photon breaks a ground state pair and creates an un- 
paired atom in state |3), (Fig. HJa)). Energy balance 
yields uj — Ej + £j. The first term is nonzero when u> 
exceeds the threshold energy ujj,. In the ground state 
6{u>) = 7T and only this term remains. The second term 
derives from excited pairs and corresponds to the process 
where a photon breaks an excited pair, (Fig. [lja)). The 
energy balance now implies u> — Ej — £j , which is nega- 
tive for all j. We see that an additional peak appears at 
lu > —Lu^p, (Fig. [21a)). The maximum in absorption is 
reached at id w — A s . Its height is suppressed, since 8(to) 
can only deviate significantly from ir in an narrow window 
of width A s around the Fermi energy, where there is a 
significant density of excited pairs. Finally, we note that 
when A s — > the two peaks merge at zero frequency, i.e. 
the RF spectrum of the gapless steady state is reversed 
to that of a normal state. 

Now let us turn to the regime of periodically oscillating 
order parameter, which occurs when (3 > tt/2 QUE!- In 
this state A(t) is given by the Jacobi elliptic function 
dn with an amplitude comparable to and a period 
of order 2n/Af. We are to analytically determine the 
Bogoliubov amplitudes using the exact solution for the 



BCS dynamics 12, yielding 
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where f2 is the frequency of oscillations of A(t), Uj — 
v(£j) is a function of single-particle energy, and 0j has 
been discussed below eq. (|4]). For brevity, the analytic 
expressions for f2, a,j n , bj n , and Vj are omitted. We note 
however that Vj plays analogous role of excitation energy 



£j for the periodic regime. For example, 
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we approach the regime of constant steady state gap, 
(3 — > 7r/2. One can also show that v(e) is a monotonic 
function of |e|, v(e) > \s\, v(0) — 0/2, and v(e) — > |e| for 
Xi — » and for large |e|. 

Comparison of the steady state wave functions ^ and 
(|4j) suggests the two terms in eq. © may be interpreted 
as two orthogonally paired states for each level j. These 
are the analog of ground state and excited pairs and have 
energies ±fj. In addition, these states contain n quanta 
of the oscillating pairing field A(t) each carrying energy 
fL We will refer to these quanta as "deltons" . These are 
quanta of the amplitude mode of the pairing field and 
can be interpreted as Higgs bosons p^l^. 

Equations (|3I6[) determine the RF spectrum in the 
periodic regime 
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FIG. 4: RF-absorption spectra eq. (J7J) for the state with pe- 
riodic in time order parameter A(t) produced by a sudden 
switch of the pairing strength A; — > Xf for Sfi ~ — 0.75A S . 
The detuning frequency ui is in units of time-averaged or- 
der parameter A s . Note two sequences of peaks at even and 
odd multiples of A a and also jumps at uj — (2n + l)A a that 
sometimes are on top of the odd peaks. The frequency of os- 
cillations of A(t) is O = 2A S . Multiple peaks are due to pro- 
cesses where an RF photon breaks an excited/ground state 
Cooper pair and emits or absorbs several oscillation quanta 
("deltons"). 



Here we dropped oscillatory terms assuming they aver- 
age to zero on the time scale of the measurement. Ex- 
pression ([7|) describes two series of equidistant peaks, 
(Fig. [4|), corresponding to the processes where an RF 
photon breaks one of the two paired states on level j and 
emits or absorbs n deltons, (Fig. [TJb)). The energy bal- 
ance reads u> — Ej ± Vj + nfi. The first series of peaks 
is described by the first term in eq. ([7|) and is analogous 
to the ground state pair peak in eq. ([5|). In this case, 
the n = peak is located at the minimum detuning fre- 
quency oj^ — v(8n) + <5/i, cf. eq. ([5]). Thus, peaks in 
the first sequence are at uj = uj^ + nfi. When A(t) is 
the Jacobi elliptic function dn, £1 = 2A S |24J, where A s 
is the time average of A(t) over the period. We note 
that the Fourier components of Bogoliubov amplitudes, 
cijn and bj n in eq. ((7]), are discontinuous at the Fermi 
level £j = similar to T = Fermi distribution. This is 
a consequence of the fact that initial states for the pe- 
riodic regime are close to the normal state H [Uj]. 
The discontinuities lead to jumps in the RF spectra at 
lu = (2n + l)A s , (Fig.0J. 

The second series of peaks is the analog of the excited 
pair peak in eq. ([5]). These peaks are at ui « — A s + nfi = 
(2n — 1)A S , (Fig. [4j. Their heights are suppressed for 
the same reason as in the excited pair peak in eq. ([5]). 
Their width is determined by the width of the excited 
pair distribution function and becomes extremely narrow 
for large (3, (Fig. [3]). In this limit, they can be superim- 
posed by jumps in the first sequence of peaks, (Fig. 2]). 

The sharp features of the RF spectra detailed above 



will be broadened by variety of effects in practice such as 
changes in particle number between experiments. More 
significant deviations will occur as one gets closer to the 
Feshbach resonance as our treatment is based on BCS 
theory. Finally, RF probing should be performed on a 
timescale shorter than the quasiparticle relaxation time 
t £ ~ Epj ' A 2 j [TBI, which limits the lifetime of the steady 
states considered here. At times larger than r e an iso- 
lated system is expected to re-thermalize to a state with 
a nonzero effective temperature which can be determined 
by balancing the total internal energy T^, 

In conclusion, we have obtained RF spectra for the 
nonequilibrium steady states formed in a fermionic con- 
densate due to a rapid switching of the pairing strength. 
The RF spectrum of the steady state with constant or- 
der parameter A s =/= has two peaks in contrast to the 
spectrum of the paired ground state where there is a sin- 
gle peak. The peak at negative detuning frequencies re- 
flects the presence of excited pairs - elementary excita- 
tions of the condensate and its shape is a direct measure 
of their distribution function. The other peak is a coun- 
terpart of the ground state spectroscopic response. In 
the steady state characterized by a periodically oscillat- 
ing A(i), each of the two peaks splits into a sequence of 
equidistant peaks with the spacing between peaks given 
by the frequency of oscillations f2. 
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